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Quantum annealing is a generic algorithm using quantum-mechanical fluctuations to search for 
the solution of an optimization problem. The present paper first reviews the fundamentals of 
quantum annealing and then reports on preliminary results for an alternative method. The review 
part includes the relationship of quantum annealing with classical simulated annealing. We next 
propose a novel quantum algorithm which might be available for hard optimization problems by 
using a classical-quantum mapping as well as the Jarzynski equality introduced in nonequilibrium 
statistical physics. 
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I. INTRODUCTION 

Recent remarkable advances in experimental tech- 
niques have enabled us to fabricate dedicated nanostruc- 
tures, in which prominent quantum effects are observed 
and/or controlled. Such technological achievements open 
up the possibilities of realization of advanced methods 
of quantum computation. In this article, we review and 
present new developments in quantum annealing (QA)i~— 
as an important subfield of quantum computation. 

Quantum annealing is a generic algorithm intended for 
solving optimization problems by use of quantum tunnel- 
ing. An optimization problem is formulated as a task to 
minimize a real single- valued function, called a cost func- 
tion, of multivariables. We will in particular be interested 
in combinatorial optimization problems, in which vari- 
ables assume discrete values^^. We will use the simple 
term of "an optimization problem" to mean a combina- 
torial optimization problem. 

Optimization problems are classified roughly into two 
types, easy and hard ones. They are distinguished by 
the computational cost to solve optimization problems. 
Easy problems are those which can be solved by the best 
algorithms in time polynomial in the problem size (poly- 
nomial complexity). For hard problems, in contrast, the 
best known algorithms cost exponentially long time as a 
function of the system size (exponential complexity) . For 
these latter problems, it is virtually impossible to find the 
optimal solution in a reasonable time, if the problem size 
exceeds a moderate value. Most of the interesting op- 
timization problems belong to the latter hard class. It 
is thus important practically to devise algorithms which 
give approximate but accurate solutions efficiently. 

One of the generic algorithms proposed as part of such 
efforts is QA. This is a method to use quantum dynamics 
and quantum tunneling effects to let the system explore 
the phase space toward the solution of an optimization 
problem. Quantum adiabatic evolution (or quantum adi- 
abatic computation) shares essentially the same idea£. 



A related algorithm of simulated annealing (SA) is 
widely used in practice to solve optimization problems 
approximately in a reasonable tim o 11 ^ 2 . In SA, the cost 
function to be minimized is identified with the energy of 
a classical statistical-mechanical system. The system is 
then given a temperature, as an artificial control param- 
eter. By decreasing the parameter slowly from a high 
value to zero, one hopes to drive the system to the state 
with the lowest value of the energy. The idea is that 
the system is expected to stay close to thermal equilib- 
rium during the protocol in SA. If the rate of decrease 
in temperature is sufficiently slow, the system will be led 
in the end to the zero-temperature equilibrium state, the 
lowest-energy state. SA is thus regarded as an algorithm 
to make use of (non-) equilibrium statistical mechanics for 
an efficient exploration of the phase space. 

Let us turn our attention to QA, which is the main 
topic of this article. In SA, we make use of thermal, 
classical, fluctuations to introduce a stochastic search for 
the desired lowest-energy state by allowing the system to 
hop from state to state over intermediate energy barri- 
ers. In QA, by contrast, we introduce non-commutative 
operators as artificial degrees of freedom of quantum na- 
ture. They induce quantum fluctuations. We first set the 
strength of quantum fluctuations to a very large value to 
search for the global structure of the phase space, simi- 
larly to the high-temperature situation in SA. Then the 
strength is gradually decreased to finally vanish to re- 
cover the original system hopefully in the lowest-energy 
state. Quantum tunneling between different classical 
states replaces thermal hopping in SA. 

An observant reader should have noticed a number 
of similarities and differences between SA and QA. In 
both methods, we have to control the relevant parame- 
ters slowly and carefully to tune the strengths of thermal 
or quantum fluctuations as desired. The idea behind QA 
is to keep the system close to the instantaneous ground 
state of a quantum system. This aspect is analogous to 
the protocol of SA, in which one tries to keep the system 
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state in quasi-equilibrium. It is indeed possible to render 
this analogy into a more precise formulation, from which 
several important developments have been achieved. An 
example of such cases will be explained later. 

Let us move our attention to some of the develop- 
ments in nonequilibrium statistical physics. Several exact 
relations between nonequilibrium processes and equilib- 
rium states have been derived. The Jarzynski equality 
(JE) is one of such remarkable results^™. This equality 
states that we can reproduce useful information on equi- 
librium states from data obtained by the repetition of 
non-equilibrium processes. We first prepare an ensemble 
of states consisting of equilibrium distribution. We then 
drive these states by tuning some parameters following 
a predetermined schedule. In general the system during 
this process cannot be kept in equilibrium. However, the 
average with weights given by the exponentiated value of 
the work done during the nonequilibrium process coin- 
cides with the thermal average of the exponential of free 
energy difference between the final and initial states in 
equilibrium. The state is not necessarily kept in equilib- 
rium also during the protocol of SA as mentioned above. 
However, JE tells us the possibility to obtain informa- 
tion of equilibrium state from nonequilibrium processes. 
Indeed, the application of JE to SA has been studied 
by several researchers^—. This application is known to 
successfully extract information on the equilibrium state 
faster than the ordinary procedure of SA with the aid of 
the property of JE. We show a quantum version of the 
algorithm with JE in this article. 

The present paper consists of five sections. In §2, we 
review the foundation of QA. In the following section, 
we introduce the classical-quantum mapping. Here we 
understand the property of SA from the point of view of 
quantum mechanics. Then, we construct an algorithm of 
QA with JE in §4. The last section presents a summary. 



II. QUANTUM ANNEALING 

Suppose that the optimization problem we wish to 
solve has been formulated as the minimization of a cost 
function, which we regard as a Hamiltonian Hq of a clas- 
sical many-body system. As a simple example, let us 
consider the problem of database search to find a specific 
item among TV items. If the specific item to be found is 
denoted by a bracket l^o), the cost function to be mini- 
mized is simply written as 

flo = /-|*o><*o|, (1) 
where I is the iV-dimcnsional identity matrix. 



A. Algorithm of QA 

In QA, we use quantum fluctuations to search the 
ground state of Hq. To this end we introduce an op- 
erator Hi, which does not commute with Hq- The whole 



quantum system is then expressed as 



H{t)= t -H +(\- t -) .//,. 



(2) 



where r is the computation time and is an important pa- 
rameter in the following discussions. The initial Hamilto- 
nian is set as H(Q) = Hi and the final system is described 
by H(t) = Hq. It is crucial to prepare an appropriate 
Hamiltonian Hi, whose ground state is trivially found. 
For instance, we prepare Hi with the following form, 



where 



Hi 



|*o)<*o|, 



(3) 



(4) 



Its ground state is trivially given by a linear combination 
of all states with a uniform weight. 

The optimization protocol by QA is to start the al- 
gorithm from the preparation of the ground state of Hi 
as the initial state. The whole system is then driven 
by the time-dependent Hamiltonian H(t) following the 
Schrodinger equation. If we choose a very large value 
of r, the system evolves slowly. Such a slow time evolu- 
tion would let the system closely follow the instantaneous 
ground state of the total Hamiltonian H{t) if the condi- 
tion set by the adiabatic theorem is satisfied. 



B. Adiabatic evolution 

Let us formulate the above statement more precisely. 
For a sufficiently large t, the adiabatic theorem guar- 
antees that the state at time t, \ip(t)), is very close to 
the instantaneous ground state |0(£)), \(0(t)\ip(t))\ w 1, if 
the following condition is satisfied, in consideration that 
the protocol of QA starts from the initial ground state. 
We assume that the instantaneous ground state |0(t)) 
is non-degenerate. If we write the instantaneous first ex- 
cited state and its energy as |l(t)} and Ei{t), respectively, 
the quantum adiabatic theorem states that condition for 
|(0(i)|V'(i))| 2 = 1 - e 2 (e < 1) to hold is given by 



:k*)i^io(*)> 



minA(t) 2 



(5) 



where max and min are evaluated between t = and 
t = t and A(t) is the energy gap, A(i) = Ei(t) — E (t) 
with E{)(t) being the instantaneous ground-state energy. 
For our case of Eq. (0 , the numerator of Eq. ([5]) with the 
time derivative of H{t) is proportional to 1/r. Therefore 
the adiabatic condition of QA is written as 



1 
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e minA(t) 2 



(6) 
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If we choose r large enough to satisfy this condition, we 
can obtain the ground state of Hq at t = r with a high 
probability. A serious problem arises if the energy gap 
A(t) collapses as the system size increases because then 
the computation time r increases according to Eq. ([5]). 
Also the requirement of a small error probability e pushes 
t to a larger value. 

We can evaluate the computation time of the database 
search problem by QA using Eqs. (dJ-Q— . The energy 
gap is given by 

Its minimum is minA(r/2) = 1/y/N. Therefore the com- 
putation time is proportional to N. This is essentially the 
same as the simple serial search. However, if we control 
the system evolution slowly around t — r/2, where the 
gap assumes its minimum, than at other time, then the 
transition to excited states is likely to be suppressed and 
the probability of success may increase. More explicitly, 
the time dependence of the Hamiltonian is chosen to in- 
volve a monotonically increasing function f(t), instead of 
the simple t/r and 1 — tjr in Eq. @, as 

H(t) = f(t)H + {l-f(t)}H 1 . (8) 

Then we demand that the adiabatic condition of Eq. ([5]) , 
without the 'max' and 'min' operations, is satisfied at 
each time t. The result is a differential equation for 
f(t), by solving which we achieve the optimal control 
of the dynamics. It actually turns out that the total 
time evolution obtained by this method is proportional 
to y/N. This is a significant improvement over the above- 
mentioned simple result t cx N. This is equivalent to the 
computation time of the Grover algorithm 22 . This is an 
explicit, remarkable example in which quantum effects 
accelerate the convergence toward the optimal solution. 
It has been established in many other cases by numerical 
and analytical methods that QA achieves superior per- 
formance compared with SA 23 . Also known is a conver- 
gence theorem of QA 2 ^: It is guaranteed that the system 
converges to the optimal solution after an infinite time of 
evolution if an appropriate control of a system parame- 
ter is chosen. This is to be compared with a convergence 
theorem for SA^i, where the time control of the system 
parameter (temperature) is constrained to be slower than 
the corresponding one in QA. 

C. Errors in QA 

Although the convergence of QA to the optimal solu- 
tion is guaranteed in the infinite-time limit as mentioned 
above, we have to carry out the computation in a finite 
time r. It is therefore important to estimate how the 
error, the probability not to reach the optimal state, de- 
pends upon t. The adiabatic theorem provides an answer 



to this question: The error after an adiabatic evolution of 
r is generically proportional to t~ 2 in the limit of large 
r as long as the system size N is kept finite 2 ^. However, 
if the system size grows, the energy gap A(t) becomes 
small and eventually collapses to zero as N — > oo when 
a quantum phase transition takes place. Thus the study 
of quantum phase transitions is an important part of the 
analyses of QA. 

For instance, how the error in QA depends on the com- 
plexity of the optimization problems is interesting in or- 
der to understand the performance of QA. The errors in 
the final state after the process of QA can be regarded as 
imperfections or defects in the state caused by the evo- 
lution across a phase transition. The dynamics across a 
phase transition is well understood via the Kibble-Zurek 
mechanism 2 -' 27 -. This aspect of error evaluation has been 
attracting increasing attention 2 ^— . 

D. Hard optimization problems 

A hard optimization problem is often represented in 
QA in terms of a Hamiltonian H(t) with a phase tran- 
sition accompanying an exponentially small energy gap 
cx exp(— aN) (a > 0). One of the examples is found in 
the random energy model^. Consider the Hamiltonian 

H = -J2 J *ui*,-,i P «---°l P , (9) 

{i P } 

where the summation is taken over all combinations 
of p sites among N of them, and of is the z compo- 
nent of a Pauli spin operator at site i. The interaction 
Ji u i 3j ... ,i follows the Gaussian distribution with a van- 
ishing mean and variance p\/(2N p ~ 1 ), which guarantees 
a sensible thermodynamic limit for any p including the 
limit p — > oo. The random energy model is defined as the 
limit p — y oo of the Hamiltonian ([9]). This limit signif- 
icantly simplifies the problem. In particular the energy 
spectrum turns out to obey the Ga ussian distribution 
with vanishing mean and variance yj N/2. Any energy 
level takes a value drawn from this distribution indepen- 
dently of other levels. 

Let us consider to find the ground state of the ran- 
dom energy model by QA. The transverse-field operator 
is often used as quantum fluctuations for searching the 
ground state in spin glasses including the random energy 
model, 

H x = -£>?. (10) 

Quantum annealing of the random energy model has the 
ground state of H\ as the initial state, which is trivially 

i—N 

given as |^(0)) = J2{a} \ a )N^ > where o, is an eigen- 
value of of. By a perturbation theory, we can evalu- 
ate the asymptotic behavior of the minimum energy gap 
of the random energy model with a transverse field as 
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A oc cxp(— aN) for a sufficiently large system size 
Therefore the adiabatic theorem states that QA, as long 
as one demands the system to stay arbitrarily close to 
the instantaneous ground state, does not work in a time 
polynomial in system size. 

We introduce another instance from one of the famous 
combinatorial optimization problems, the exact cover. 
This is a special version of constraint satisfaction prob- 
lems. Let us consider TV spins and M "clauses" , each one 
of which involves three Ising spins chosen at random. The 
energy of a clause is zero if one spin is —1 and the other 
two are 1. Otherwise, the energy is defined as 1. This 
assignment of the energy is realized by an Ising model as 



rn—l 

+ cr m 1 cr m 2 + + ^maS + 5 J • 

(11) 

Numerical simulations have shown that the minimum en- 
ergy gap during QA of the exact cover has the asymptotic 
form proportional to exp(— aJV) 35 ' 36 . Thus we again en- 
counter a similar problem to the random energy model. 

The behavior of the energy gap as exp(— aN) is char- 
acteristic of first-order phase transitions in quantum sys- 
tems. Therefore the problem concerning the bottleneck 
of QA appears even in the simple ferromagnetic model, 
namely the Ising model with three-body interaction^. 
This model also exhibits a first-order phase transition 
with an exponentially small gap. This fact may be taken 
as an indication that QA would fail to find a trivial 
ground state in a reasonable time. However we have some 
room in the choice of quantum fluctuations expressed by 
H\. It remains to see if an appropriate choice of the 
quantum term Hi can help us avoid these difficulties in 
the framework of QA under the constraint of an adiabatic 
evolution. 

In the rest of this article, we provide a different ap- 
proach to overcome the difficulty of the energy gap clo- 
sure based on a different idea. 



III. CLASSICAL-QUANTUM MAPPING 

The analogy between QA and SA can be pursued in 
a systematic manner by the exploitation of an ingenious 
mapping between a classical thermal state and a cor- 
responding quantum stated. This mapping is differ- 
ent from the conventional quantum-classical mapping in 
terms of the Suzuki- Trotter formula or the path integral 
formulation. In particular, the same spatial dimensional- 
ity is shared by classical and quantum systems according 
to the present mapping. 



A. Mapping of a classical system to a quantum 
system 

In order to formulate the classical-quantum mapping, 
we consider a classical Ising spin system with a config- 
uration a. This system is assumed to obey a stochastic 
dynamics described by the master equation 

j t P{u;t) = Y,M{a\a';t)P{a';t), (12) 

cr' 

where M(a'\a;t) expresses the transition matrix satisfy- 
ing the detailed balance condition 

M(cr|o-'; t)P(a') = M(a'\a; t)P(a). (13) 

The instantaneous equilibrium distribution 
e -««)ff«W/2(() is denoted as P(a;t). Since we 
are considering a dynamic control of the temperature, 
T(t), the time variable t has been written explicitly 
in the arguments of the inverse temperature and the 
partition function. 

The (instantaneous) equilibrium state of the above 
stochastic dynamics is now shown to be expressed as the 
ground state of a quantum Hamiltonian. A general form 
of such a quantum Hamiltonian is given as 

H q (a'\a;t) = 5 a , a - e W H °W 2 M(a'\*;t)e-W H oW 2 . 

(14) 

This Hamiltonian has the following state as its ground 
state, 

|* cq (*)) =-i=^e-^Wk). (15) 

V Z \ t ) a 

It is clear that the quantum expectation value of a physi- 
cal quantity A(a) by |\E' oq (t)) coincides with the thermal 
expectation value of the same quantity. The ground state 
energy for this ground state is 0, which is shown by the 
detailed-balance condition, 

(^,,-e^^)M(aV;t)e-^ ff »W) |* eq (i)) 

oc £ ( e -^ - e^ H ^M(a\a'; t)) \a) 

a 

= 0, (16) 

where we used the detailed balance condition and 
the requirement for conservation of probability 
M(a\a'; t)\a) = \a"). The excited states have 
positive-definite eigenvalues, which can be confirmed 
by the application of the Perron-Frobenius theorem. 
By using the above special quantum system, we can 
treat a quasi-equilibrium stochastic process in SA as an 
adiabatic quantum-mechanical dynamics in QA. The 
stochastic dynamics can be mapped into the quantum 
dynamics as above introduced even for the system 
with continuous variables by focusing the relationship 
between the Fokker-Planck equation and the Shrodinger 
equation^. 
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B. Similarity of SA and QA 

Let us consider an optimization problem that can be 
expressed as a Hamiltonian with local interactions 



(17) 



where Hj involves cj 



and a finite number of <r 
j). A simple instance is a system with nearest- neighbor 
interactions, 



H 



(18) 



k£j 



where k is a site neighboring to j. The following explicit 
choice of the quantum Hamiltonian as an example of Eq. 
(TlH) facilitates our analysis, 



'-H l 



(19) 



The coefficient c is exponentially small in N . The above 
result reproduces the Geman-Geman condition for con- 
vergence of SASi. These authors used the theory of 
the classical time-dependent Markov chain representing 
nonequilibrium stochastic processes. The classical sys- 
tem under such a time evolution may not seem to stay 
close to equilibrium since the temperature changes con- 
tinually. However, the result mentioned above implies 
that the quasi-equilibrium condition is actually satisfied 
during the process following Eq. (j2"2")l because the adi- 
abatic condition means that the quantum system is al- 
ways very close to the instantaneous ground state, which 
represents the equilibrium state of the original classical 
system. 

The classical-quantum mapping explained in this sec- 
tion will be shown to be an important ingredient of the 
development in a later section. 



IV. JARZYNSKI EQUALITY 



where = e~^(*) p withp = maxi\Hi\, which is propor- 
tional to the energy and is of O(N ) because of finiteness 
of the interaction range. 

The quantum system (fl^)) has the following interest- 
ing property. For very high temperatures, f3 — > 0, this 
quantum Hamiltonian reduces 



(20) 



The ground state of Eq. (I2U)) is the uniform linear com- 
bination of all possible states in the basis to diagonalize 
{erf}, which means that all states appear with an equal 
probability. Such a situation is realized also in the high- 
temperature limit of the original classical system repre- 
sented by Eq. (fT7| . Thus the quantum ground state 
in the limit f3 — > appropriately describes the classical 
disordered state. Analogously, in the limit (3 — > oo, the 
quantum Hamiltonian becomes 



#nt)~x(*) 



H, 



(21) 



whose ground state is the ground state of the classical 
system (flT]) because each Hj takes its lowest value. It 
has therefore been confirmed that the limit j3(t) — > oo 
of the quantum system correctly corresponds to the low- 
temperature limit of the classical system. 

Surprisingly, the adiabatic condition applied to the 
quantum system (fll?]) can be shown to lead to the condi- 
tion of convergence of SA22,. If we assume a monotonic 
decrease of the temperature T(t) = l//3(i), the adiabatic 
condition, Eq. ([5]) without max and min, applied to the 
quantum Hamiltonian (|19l) yields the following time de- 
pendence of f3(t) in the limit of large N, 



log(cf +1) 
pN 



(22) 



We now develop a new algorithm for optimization 
problems using the classical-quantum mapping and the 
Jarzynski equality. This latter equality relates quantities 
at two different thermal equilibrium states with those of 
non-equilibrium processes connecting these two states. 
It can also be regarded as a generalization of the well- 
known inequality of thermodynamics, AF < (W). Here 
the brackets (• • ■ ) are for the average taken over non- 
equilibrium processes between the initial and final states, 
both of which are specified only macroscopically and thus 
there can be a number of microscopic realizations. 



A. Formulation of the Jarzynski equality 



The Jarzynski equality is written as 

^0 



(23) 



Here the partition functions for the initial and final 
Hamiltonians are expressed as Zq and Z T , respectively. 
One of the important features is that JE holds indepen- 
dently of the pre-determined schedule of the nonequi- 
librium process. Another celebrated benefit is that JE 
reproduces the second law of thermodynamics given as 
AF < (W) by using the Jensen inequality, {e~^ w ) > 
e -P(w) _ Notice that we have to take all fluctuations into 
account in formation of the expectation value in the right- 
hand side of JE, if one uses the calculation of the free 
energy difference in practical use. Finite errors on the 
ratio of the partition functions are attributed to the rare 
event on nonequilibrium processes^—. 

We show the proof of JE for a classical system in 
contact with a heat batbM. Let us consider a thermal 
nonequilibrium process in a finite-time schedule < t < 
t. Thermal fluctuations are simulated by the master 
equation. We employ discrete time expressions and write 
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tk+i — t k — St, to — and t n = r. The probability that 
the system is in a state 07. at time tk will be denoted 
as P{(Tk\tk)- Notice that a k is not a state of a single 
spin at site k but is a collection of spin states at time 
tk- The transition probability per unit time is written 
as M(ak+i \c k ', tk)- In the original formulation of JE, the 
work is defined as the energy difference due solely to the 
change of the Hamiltonian. We can construct JE in the 
case of changing the inverse temperature, which is use- 
ful in the following application, by defining the work as 
-f3W(a k ;t k ) = -W(t k+1 ) - P(t k ))E(a k ), where E{a) 
is the value of the cost function (classical Hamiltonian 
Hq((j)) with the spin configuration a. The left-hand side 
of JE is then expressed as 



zn{ 

{a k } fe=0 

xP(<7 ;t ) 



e -t3W{a k+1 ;t k ) e StM(a k + 1 \(j k ;t k ) 



(24) 



This means that the system evolves from state a k to state 
CTfc+i during the interval St and then we measure the work 
W(a k+ i;t k ). This process is repeated from k = to 
n — 1. The initial condition is set to the equilibrium dis- 
tribution. If the transition term exp(StM(a k+ i\cr k ;t k )) 
is removed in this equation, JE is trivially satisfied be- 
cause the summation of —/3W(a k +i;tk) over k yields 
— (/3(i„) — (3(to))E(ao)). A non-trivial aspect of JE is 
that the insertion of the transition term does not alter 
the conclusion. 

Let us evaluate the first contribution k = of the 
product in the above equation as 



l e -PW(<7i;to) e 6tM (<Ti\<7 ;to) 



}P{<ro;t ) 



= _ e -G8(ti)-0(t o ))£(<ri)g-0(ta)E(ffi) 



Zi 



(25) 



Repeating the same manipulations, we obtain the quan- 
tity on the right-hand side of Eq. (|2"4"|) as 



71-1 _ 

^k 



k=0 



Zq 



(26) 



exp(-(/9(i k+1 ) - f3(tk))H (a k )/2). If we apply this oper- 
ation to the quantum wave function |5' cq (to)), the state 
is changed into a state corresponding to the equilibrium 
distribution with the inverse temperature /3(t k +i). Even 
if the quantum time-evolution operator U(a'\<j) tk+i) — 
exp(—iH q (a'\a;tk+i)St) is applied, this state does not 
change, since it is the eigenstatc of H q (a'\a;tk+i)- The 
obtained state after the repetition of the above procedure 



is given as 



n-l 



k=0 



|*(*n)) oc ~[[{W(a k+1 ;t k )U k+1 (a k+1 \a k ;t k )} 

Ee-^IcA (27) 



B. Quantum Jarzynski Annealing 



This is essentially of the same form as Eq. (|24p . In- 
stead of M(a k+ i\cr k ;t k ), we use the time-evolution oper- 
ator U(cr k+1 \a k ;t k ) here. After the system reaches the 
state \^(t n )), we measure the obtained state by the pro- 
jection onto a spin configuration a'. The probability is 
then given by |(cr'|\E'(t rl ))| 2 . The ground state is obtained 
by the probability proportional to exp(— f3(r)H ), since 
|\&(f n )) oc I 'I'cq (£«))■ If we repeat the above procedure 
for (3(t) 3> 1, we can efficiently obtain the ground state 
of Hq. This is called the quantum Jarzynski annealing 
(QJA) in this paper. 

It may seem unnecessary to apply the time-evolution 
operator U(a k +i\cr k ; t k ) at the middle step between the 
operations of work exponentiated operators W(a k +i; t k ). 
In the formulation of JE, we can use the identity matrix 
instead of U(a k +i\cr k ; t k ), but this is not perfectly con- 
trollable in quantum computation. The quantum system 
always changes following the Shrodinger dynamics. The 
change by quantum fluctuations during the protocol is 
described by the time evolution operator U(a k +i \<J k ; t k ) 
at the middle step between the operations of work expo- 
nentiated operators W{<Jk+i] tk)- However, it is useful to 
remember the nontrivial point of JE. When transitions 
by quantum fluctuations between the exponentiated work 
occurs, JE holds as in Eq. (|2"3|). 

We emphasize that the scheme of QJA does not rely 
on a quantum adiabatic control. Aside from the number 
of the unitary gates to build the quantum circuit of QJA, 
the computation time does not depend on the energy gap. 
Therefore QJA would not suffer from the energy-gap clo- 
sure differently from ordinary QA. The result should be 
independent of the schedule of tuning the parameter j3(t) 
in the above manipulations. 



Let us consider the application of the above formula- 
tion of JE to QA using the quantum system (|f4[) . Ini- 
tially we prepare the trivial ground state with a uniform 
linear combination as in ordinary QA. From the point 
of view of the classical-quantum mapping, this initial 
state expresses the high-temperature equilibrium state 

|*eq(*o)) OC e -P(to)H a (a)/2^ ^ < J Then w& 

introduce the exponentiated work operator W(a k ;tk) = 



C. Performance of QJA 

We expect that QJA gives the probability of the 
ground state following the Gibbs-Boltzmann factor un- 
der any annealing schedule. In contrast, without the 
multiplication of the exponentiated work, slow quantum 
sweep is necessary for an efficient achievement to find 
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the ground state according to the ordinary QA. Let us 
consider two instances of the application of QJA. 

The first example is the database search to find a spe- 
cific item among N = 50 items as explained in Sec. II. 
We reformulate this database search as a single mini- 
mum potential among uniform values as Vq = —1, and 
Vi = (i 0). We construct the transition matrix fol- 
lowing the Glauber dynamics. We set up QA by the 
classical-quantum mapping and QJA by use the same 
transition matrix for comparison. We increased j3(t) from 
to 10 linearly as a function of t. Figure Q] describes the 
results by QA and QJA for different computation time 
t = 1,2 and 5. The plots by QJA (upper curves) are 
indistinguishable from the reference curve representing 
the instantaneous Gibbs-Boltzmann factor for the mini- 
mum state i — 0. In contrast, one finds that QA (lower 
curves) needs sufficiently slow control of the quantum 
fluctuations to efficiently find the minimum. 

The second example is to search the minimum from a 
one-dimensional random potential. We consider to find 
the minimum of the Hamiltonian H$ = — X}j=i 
where Vi denotes the potential energy at site i. We choose 
{Vi} randomly. We again compare the case for TV = 50 
sites by QJA and QA. By a linear schedule for increasing 
the parameter j3(t) from to 100, we carry out QA with- 
out the exponentiated work operations and QJA. Figure 
[2] shows the comparison between the probability for find- 
ing the ground state by QA and QJA for different sched- 
ules t — 1, 10, and 100. On finds that the plots by QJA 
(upper curves) are again indistinguishable from the ref- 
erence curve. QA (lower curves) needs sufficiently slow 
decrease of quantum fluctuations to efficiently find the 
site with the minimum potential. 

QJA does not suffer from the energy-gap closure by in- 
creasing the system size and gives the probability to ob- 
tain the ground state independently of a predetermined 
schedule. Both of the above cases indeed show that the 
results of QJA does not depend on the computational 
time t. 

To realize QJA in an actual quantum computation, 
we need to implement the exponentiated work operation 
W(<Tk',tk), since it looks like a non-unitary operator. We 
can construct this operation as a unitary operation by us- 
ing an extra quantum syste m 43 ' 44 . Instead of cut "time" , 
we need another resource like a "memory" in quantum 
computation. Fortunately, the amount of needed mem- 
ory does not diverge exponentially by an increase of the 
system size N. Thus the results by QJA shown here 
imply that we may overcome the difficulties in hard opti- 
mization problems and solve them in a reasonable time. 
However, the above instances do not represent hard op- 
timization problems as in the preceding section with an 
exponentially vanishing energy gap. In that sense, our re- 
sults are preliminary and future work should clarify the 
efficiency of QJA for such hard problems. 
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FIG. 1: Quantum Jarzynski annealing for the database search 
problem from N = 50 items (N denotes the size of the Hilbert 
space) . The probabilities to find the desired state are plotted 
for r = 1, 2, and 5 from top to bottom. The dashed curves 
denote the instantaneous Gibbs-Boltzmann factor for refer- 
ence. The upper solid curves (blue curves) representing the 
results by QJA are fixed to these reference curves, while the 
lower ones (red ones) express the time-dpendent results by 
the ordinary QA. 



V. SUMMARY 

In this paper we first reviewed QA, which is a generic 
algorithm to approximately find the solution to optimiza- 
tion problems. The method is based on the adiabatic 
evolution of quantum systems. A sufficient condition to 
efficiently find the ground state by QA is given by the 
inverse square of the energy gap. The same analysis in 
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conjunction with the classical-quantum mapping enables 
us to reproduce the convergence condition of SA, which 
is another approximate solver exploiting thermal fluctu- 
ations. Recent investigations disclosed the bottleneck of 
QA, the exponential closure of the energy gap. Accord- 
ing to the prediction by the adiabatic theorem, QA does 
not work efficiently on problems with such a bottleneck. 
We thus have to invent an alternative way to efficiently 
find the optimal solution with the energy gap vanishing 
exponentially with the system size. 



We next showed an application of JE to QA as one 
of the possible improvements of QA, using the fact that 
we can map a quantum system close to the instanta- 
neous ground state to a thermal classical system in quasi- 
equilibrium by the classical-quantum mapping. This pro- 
tocol keeps the quantum system to be the ground state 
expressing the instantaneous equilibrium state. The cost 
for the realization of QJA in quantum computation does 
not diverge exponentially as the system size increases, 
which is the essentially different point from ordinary QA. 
This idea has been shown to work in a few examples with 
a small Hilbert space. It is a future problem to test this 
method for problems of much larger scales and hard ones 
with the exponential closure of the energy gap. 
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FIG. 2: Quantum Jarzynski annealing for the random- 
potential problem for N = 50 sites (N denotes the size of 
the Hilbert space). The probabilities to find the minimum- 
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bottom. The same symbols as in Fig. fjjare used. 
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